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We present the nonlinear susceptibility as a direct test of the quadrupolar Kondo 
scenario for heavy fermion behavior, and apply it to the case of cubic crystal-field 
symmetry. Within a single-ion model we compute the nonlinear susceptibility result- 
ing from low-lying (5/ 2 ) and Kramers (5/ 3 ) doublets. We find that nonlinear sus- 
ceptibility measurements on single-crystal U Be±3 are inconsistent with a quadrupolar 
(5/ 2 ) ground-state of the uranium ion; the experimental data indicate that the low- 
lying magnetic excitations of UBe\z are predominantly dipolar in character. 
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There exist several metallic systems whose novel thermodynamic, magnetic and trans- 
port properties are not adequately described by conventional Fermi liquid theory; specific 
examples include the quasi-one dimensional conductors, JIJ certain actinide heavy fermion 
materials 0-0 and the layered cuprate superconductors. [J7||| The search and characteri- 
zation of non-Fermi liquid (NFL) fixed points is thus a topic of active research. [^-0 In 
this Letter we present an unambiguous experimental test of the quadrupolar Kondo effect, 
a model proposed by Cox || to characterize the NFL behavior observed in the cubic three- 
dimensional heavy fermion material UBei^. We use the nonlinear susceptibility (^3) as a 
direct probe of low-lying quadrupolar fluctuations, and compute its behavior within a single- 
ion model for the case of cubic crystal-field symmetry; these predictions are then compared 
to X3 measurements on single-crystal UBe^. 

Most heavy fermion metals display a dramatic reduction in resistivity at low temperatures 
associated with the development of coherent quasiparticle propagation. UBe^ is atypical, 
undergoing a superconducting transition directly from a normal state with a large incoherent 
resitivity || of order 140/if2 cm. The low-temperature dependences of the magnetic suscep- 
tibility and the specific heat are logarithmic in the approach to the superconducting 
transition. Resistance, |]TU| specific heat, [[□]], susceptibility [0 and magnetoresistance |T3| 



measurements indicate that Fermi liquid behavior is restored at low temperatures under an 
applied pressure. Thus UBe^ is a metal with a tunable Fermi temperature (Tp) such that 
Tp < T c at ambient pressure. The microscopic physics underlying this suppressed, pressure- 
dependent [pT0| p~3j| and field-dependent |14| Tp is a crucial issue for the characterization of 



the complex many-body ground-state of UBe^. 

Cox |jj has proposed that novel single-ion physics is responsible for the observed NFL 
behavior in UBe^. The observation of a well-defined Schottky anomaly |15[ at T ~ 180K 
indicates that the uranium ion is in a local-moment rather than an intermediate valence 
regime; |16| however, it can assume either the U 4+ (5/ 2 ) or the U 3+ (5/ 3 ) nominal valence 
state, [|5|,|T51| and neither quasi-elastic neutron scattering JT7J nor photoemission ]T8| , |TP 



mea- 



surements can unambiguously resolve the crystal-field assignments. In a cubic enviroment 



Cox has identified a non-magnetic quadrupolar (T3) ground-state p|,pO| of the U 4+ ion. He 
suggests that fluctuations within this non-Kramers doublet are overscreened by the con- 
duction electrons; this quadrupolar Kondo effect then leads naturally to a non-Fermi liquid 



ground-state. JzlJ] In an alternate scenario, supported by NMR measurements JsIJ consistent 
with a U 4+ valence state, the low-lying spin excitations are dipolar; the NFL behavior is 
attributed to the system's proximity to a T = quantum phase transition, [|]|| analogous 
to that recently observed in MnSi by Lonzarich and coworkers. [^3] 

The nonlinear susceptibility (^3) is an ideal test of the quadrupolar Kondo effect in 
UBeis; it can distinguish unambiguously between a low- lying quadrupolar and Kramers 
crystal-field doublet. In the paramagnetic state, %3 measures the leading nonlinearity in the 
magnetization 



M = X iB + ^XsB 3 + ... 



(1) 



in the direction of the applied field (B); it was originally proposed as a direct probe of order- 
parameter fluctuations in spin glasses. []23J Morin and Schmitt |25[ extended this technique to 
non-random spin systems, where they used the nonlinear susceptibility to study quadrupolar 
interactions in rare-earth intermetallic compounds. 



The most general form for xs m a cubic enviroment is |26| 



Xs = Xs + AX3$(&) 



(2) 



where $(&) is the cubic harmonic 



$(&) = - 3(6 



K + K 



1 



(3) 



and the bi (i = 1, 2, 3) are the direction cosines of the field. The numerical factors in $(6) 
are chosen so that A% 3 = xl 00 ~ Xz j the "powder-averaged" component of the nonlinear 
susceptibility is X3 = xl U + ^ A X3- 

The ratio of the two contributions to X3 m (2) is qualitatively different for a quadrupo- 
lar and a magnetic ground-state. An isolated Kramers doublet |25j results in a nonlinear 



susceptibility X3 — ~ir5 that is isotropic, reflecting the negative curvature of the Brillouin 



function. For an isolated doublet with a quadrupolar moment Q, the field-dependent part 
of the Hamiltonian is H = \B 2 Q ab b a b bl where Q a b oc [J a Jb — \5 a bJ(J + 1)] is the quadrupole 
operator; [EH] more explicitly 

QB 2 \ I** Qx* ~ 
Qxx + kyy ~Q 



H 



1 z z 



(4) 
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where q aa = b 2 — | (a — x,y, z). Diagonalizing H, we find that the splitting of the quadrupo 
lar doublet is given by 



El ± = E T ±QB 2 ^ (5) 

which yields an anisotropic nonlinear susceptibility Xz(b) = zf*&(b)- In Cox's model for 
UBeis, there is partial quenching of Q by the conduction sea and 



Q 2 J % T»T n 

21 ' ^rMTo/T) T«T 



where Tq is the "Bethe ansatz" Kondo temperature; [^7| the exact solution of the two 



channel Kondo model [p7|-p9| yields an asymptotic form for f(x) with an associated value 
f2?|j a = 1/tt 2 w 0.10. Thus the quadrupolar Kondo hypothesis predicts a Ax3{b) that 
increases logarithmically with decreasing temperature. 

In this idealized discussion we have neglected the Van Vleck contributions to X3- I n 
practice, a uranium atom in a magnetic configuration {U 3+ ) with a moment fi(H) = fiQ + 
B will exhibit a small A% 3 = 4 ^°^( fe ) cl ue to the nonlinearity in fi(H). Conversely, 
X3 for a U ion with a low-lying quadrupolar doublet {U A+ ) will have an isotropic Van 
Vleck component {xX V ) that has a weak temperature dependence. Despite these additional 
contributions, we expect the anisotropic component of the nonlinear susceptibility 

Ax 3 (fe) = (x 3 (fe) - *n m 

to be small and nearly temperature-independent for a uranium atom with a dipolar ground- 



A X3 (fc) 

X'i 

fluctuations are quadrupolar in nature. 



state; by contrast, should be large and strongly temperature-dependent if the low-lying 



Single-ion crystal-field calculations allow us to quantify the preceeding discussion, and 
they have been performed for J = 4 and J — | manifolds of / orbits in a cubic enviroment. 
PDfl The overall energy scale (W) and the level ordering (x) have been adjusted to fit the 
total entropy in the observed Schottky anomaly, ]nj and the resulting energy schemes are 
displayed in Figure 1. The associated Xz(b) and Axs ^ are shown in Figure 2 and 3 respec- 

X3 

tively, where the moment has been normalized by a fit to the measured high-temperature 
susceptibility; |HJ the numerical solution of the two-channel Kondo model |27j has been 



used to determine the effects of screening in Figure 3. 

In order to test the quadrupolar scenario in UBei3, we measured X3 along the three 
principal crystal axes of an oriented single crystal grown from Al flux. The superconducting 
transition temperature, a rough measure of the sample quality, was found by specific heat 
to be T c = 0.75K for this crystal. Measurements were also performed on a a polycrystalline 
sample with T c = 0.96. Finally a third sample, an unoriented single crystal with T c = 0A8K, 
was studied. For the X3 measurements on the oriented crystal, the orientation was achieved 
with a precision of ±3 degrees; the data were taken as M vs. B at fixed temperatures up 
to 4 Tesla in a Quantum Design SQUID magnetometer. The deviation from linearity was 
only ~ 2% at the lowest temperature and the highest field; it was attributed to the leading 
nonlinear contribution of M to ^3- The magnetization data were fit to the expression 
M = Xo + XiB + ^X3B 3 , where Xo was included to avoid systematic errors associated with 
both trapped flux in the superconducting solenoid and a small (~ 10 ppm), ubiquitous 
ferromagnetic signal which saturated at ~ 1 Tesla. The temperature dependence of X3 is 
displayed in Figure 4. The data were typically fit over the region 2 < B < 4 Tesla; in this 
field range, M/B — Xo wa s always linear with respect to B 2 . The linear part, Xi ( n °t shown), 
agrees well with published values. || Figure 4 shows the nonlinear susceptibility measured in 
the 111, the 110, and the 100 directions. We note that the observed X3 ls both negative and 
monotically decreasing with decreasing temperature; its magnitude is significantly greater 
than that predicted for the quadrupolar scenario (Figure 2a), but comparable in size at 
T ~ 1§K with that expected from a dense concentration of only partially quenched U 



magnetic doublets. 

The observed magnitude and temperature dependence of X3 was similar for the other 
two samples studied. The measurements on the polycrystalline sample (Figure 4) provide a 
crucial control on our results; here we expect the impurity level to be low given the relatively 
high observed value of T c . The polycrystalline sample displays behavior in similar 
to that of the orientation-averaged single- crystal. This result, combined with the large 
magnitude of X3, exclude the possibility that the observed %3 is due a residual background 
of magnetic impurities. 

The measured anisotropy in the nonlinear susceptibility (Figure 4 inset) is small ( Ax3 ^ 6 ) ~ 
3 x 10 _1 ) with a very weak temperature dependence; moreover it appears, at the level of one 
standard deviation, to have the opposite sign to that expected for the quadrupolar scenario 
(see Figure 3). These results strongly favor a magnetic model for the uranium ions in UBe^ 
with a low Kondo temperature. One can try to reconcile these results with the quadrupolar 
scenario by invoking a large Van Vleck contribution (x^ V )] it would result from virtual spin 
or valence fluctuations into higher lying multiplets of the U ion. Such a term would scale 
approximately with where A x is the gap to the higher multiplets. In order for x\ v ~ ^- 
to be much larger than A% 3 ~ ^ we need < T , a condition inconsistent with the initial 
assumption of a well-defined quadrupolar ground-state. 

We now return to the possible origins of NFL behavior in UBe^. Though a single-ion 
mechanism cannot be ruled out, j31|a canonical Kondo model for the magnetic U ion results 



in a Fermi liquid ground-state. Furthermore one expects a system with a low-lying Kramers 
doublet to display a reduction in 7 = ^ when g^isB ~ T F , in contrast to that observed 
|32|1 for UBe±3. Thus we conclude that these results cannot be explained within a singe-ion 



picture picture and require a more sophisticated approach, possibly one that has an intrinsic 
pressure- and field-dependent T F . We are tempted to identify the observed NFL behavior 
as a lattice phenomenon, possibly attributed to the system's proximity to a T = quantum 
phase transition. Two different types of experiments would clarify this situation. First, 
thermodynamic and transport studies on U x Th\- x Be\z would probe the behavior of dilute 
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U atoms in the cubic enviroments, |33| thereby indicating the importance of lattice effects. 
Second the nonlinear susceptibility as a function of pressure could be measured; we expect a 
shoulder in X3 ~ w that coincides with the observed development of Fermi liquid behavior 
in the resistance, [TIJ the specific heat [D], and the magnetoresistance. JT3[ 

In conclusion we have performed a series of nonlinear susceptibility measurements on 
the cubic heavy fermion system UBe^. We find a small weakly temperature-dependent 
anisotropy, ~|~ , in the nonlinear susceptibility that is difficult to reconcile with the 
quadrupolar Kondo scenario. These results provide strong evidence for a Kramers doublet 
ground-state in the U 3+ ions of UBeis and suggest a lattice mechanism for the observed 
non-Fermi liquid behavior. Further experiments have been proposed to test this conjecture. 

We thank D.L. Cox and A.M. Tsvelik for extensive discussions related to this work. P. 
Coleman is supported by NSF grant DMR-93-12138 and work at Los Alamos was performed 
under the auspices of the UCDOE. 
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Fig. 1. The J = 4 quadrupolar and the (b) J = | dipolar single- ion energy schemes 
for UBeiz where the overall energy scale and the level ordering are determined by a two- 



parameter fit to the specific heat measurements of Felten et al. [15 
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Fig. 3. The anisotropic part of the nonlinear susceptibility for the J = 4 level 
scheme of Fig. 1. The low-temperature (dotted line) was determined by 

X3 

normalizing the single-ion anisotropy with the screening function /(T/Tq) from 



the solution of the two-channel Kondo problem; ||27|| here the value T = 1.5K 



was extracted from the observed specific heat [|i~5[1 . For the J = | scheme of 



Figure 1 
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Fig. 4. The measured nonlinear susceptibility (x3(b)) and Ax3< ^ (inset) for single 
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crystal and polycrystalline UBe lz . 
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